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KEY  PHYSICS  ISSUES  AFFECTING  THE  PERFORMANCE 
OF  FREE  ELECTRON  LASERS 

1  INTRODUCTION 


Free-electron  lasers  (FELs)  comprise  a  class  of  potentially  efficient  devices 
capable  of  generating  high  quality  coherent  radiation,  continuously  tunable  from 
the  microwave  region  through  to  the  ultraviolet  region,  at  high  average  and/or 
high  peak  power.1-20  Although  the  first  ideas  for  the  FEL  have  been  known 
since  the  original  article  by  Motz21  in  1951  and  the  successful  experiments  by 
Phillips22  in  1960,  it  is  the  work  by  Madey23-25  in  the  infrared  wavelength 
regime  in  the  1970s  that  has  made  the  FEL  a  serious  candidate  for  a  new  pow¬ 
erful  radiation  source. 

The  FEL  is  characterized  by  a  relativistic  electron  beam  co-propagating 
through  an  undulator  field  with  an  input  radiation  field.  If  the  electron  beam 
is  of  good  quality  and  has  a  sufficiently  high  current  density,  input  radiation 
at  the  appropriate  frequency  will  grow,  as  in  the  case  in  an  FEL  amplifier. 
The  undulator  field  can  be  a  static  periodic  magnetic  field,  a  static  periodic 
electric  field,  or  an  electromagnetic  radiation  field.  When  the  gain  per  pass  is 
low  through  the  undulator,  many  passes  are  required  for  the  radiation  to  reach 
saturation,  as  in  an  FEL  oscillator.  A  schematic  of  the  FEL  oscillator  is  shown 
in  Fig.  1.  Here  the  end  mirrors  reflect  the  radiation  between  the  exit  and  the 
entrance  of  the  undulator. 

The  FEL  has  major  differences  relative  to  conventional  lasers.  The  active 
medium  in  the  FEL  is  not  atoms  or  molecules  but  free  electrons.  The  radiation 
frequency  of  the  FEL  is  continuously  tunable  and  i6  not  limited  to  the  discrete 
allowed  quantum  transitions.  The  free-electron  lasing  medium  also  has  the  ad¬ 
vantage  that  it  can  support  large  amplitude  electromagnetic  fields  that  could 
damage  conventional  material  lasers.  In  the  FEL,  the  electron  kinetic  energy  is 
converted  into  radiation.  Some  electrons  lose  kinetic  energy,  while  others  gain 
kinetic  energy.  Under  the  appropriate  conditions,  coherent  radiation  is  gener¬ 
ated  as  electrons  on  the  average  lose  kinetic  energy.  Since  each  electron  can 
gain  or  lose  many  photons,  the  physics  can  be  understood  in  terms  of  classical 
mechanics. 

Electrons  can  lose  kinetic  energy,  as  they  travel  through  the  interaction  re¬ 
gion,  in  two  ways,  i)  Due  to  electron  acceleration  by  the  undulator  field  alone, 
magnetic  bremsstrahlung  radiation  is  produced.  When  the  the  electrons  are 
initially  uniformly  distributed  axially  in  the  beam  (not  bunched  with  a  peri¬ 
odic  wavelength),  this  spontaneous  radiation  is  incoherent  and  low  in  power, 
ii)  The  electrons  can  lose  significant  kinetic  energy  due  to  the  force  from  the 
“ponderomotive  potential  well,”  which  is  a  traveling  beat  wave  formed  by  the 
radiation  and  the  undulator  fields.  When  the  axial  electron  beam  velocity  is  ap¬ 
proximately  equal  to  the  phase  velocity  of  the  ponderomotive  potential  well,  the 
initially  uniformly  distributed  electrons  become  trapped  in  the  ponderomotive 
potential  well  and  bunched  at  the  ponderomotive  wavelength.  The  stimulated 
radiation  from  the  periodic  beam  bunches  is  coherent  and  is  many  orders  of 
magnitude  higher  in  power  than  the  spontaneous  radiation. 
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The  key  physics  issues  concerning  the  generation  of  the  stimulated  coherent 
radiation  are  the  focus  of  this  paper.  These  issues  include,  among  others, 

i)  generation  and  propagation  of  high  quality,  high  current  electron  beams, 

ii)  optical  guiding  of  the  generated  radiation  beam, 

iii)  undulator  field  errors,  and 

iv)  excitation  of  sideband  radiation. 

Electron  beam  sources  vary  widely  in  energy,  beam  quality,  current,  pulse 
length,  and  repitition  rate.  Radio  frequency  (rf)  driven  accelerators,  such  as  rf 
linacs,  microtrons,  and  storage  rings,  have  peak  currents  from  a  few  amperes 
to  hundreds  of  amperes,  electron  energies  up  to  hundreds  of  MeV,  and  electron 
micro-pulse  lengths  much  shorter  than  the  rf  wavelength.  Induction  accelerators 
can  have  many  kilo-amperes  of  current  and,  currently,  have  energies  up  to  tens 
of  MeV  and  pulse  lengths  of  a  few  tens  of  nanoseconds  to  a  few  microseconds. 
Electrostatic  accelerators  usually  have  low  current  and  low  energy,  but  good 
beam  quality.  Lumped  parameter  modulators  usually  can  provide  relatively 
inexpensive  electron  beams  in  the  low  energy  range  with  pulse  durations  up  to 
several  microseconds  and  hundreds  of  amperes  of  current.  Electron  beams  from 
pulsed  power  sources  can  have  very  high  currents,  but  they  have  low  energies 
and  pulse  durations  of  tens  of  nanoseconds.  Each  of  these  electron  beam  sources 
have  been  applied  to  the  generation  of  FEL  radiation.  The  output  radiation 
characteristics  vary  widely  as  a  result  of  this  wide  variation  of  electron  beam 
sources. 

Free-electron  laser  experiments  are  conducted  all  over  the  world:  USA, 
France,  Italy,  Germany,  USSR,  United  Kingdom,  the  Netherlands,  China,  Japan, 
and  Israel.  A  recent  compilation  of  the  FEL  projects  around  the  world  is  given 
in  Ref.  20.  A  brief  summary  of  FEL  facilities  in  the  United  States  is  given  here. 
Although  the  conceptual  configuration  of  the  FEL  is  simple,  the  implimenta- 
tion  was  arduous  because  existing  electron  beam  sources  were  not  designed  with 
FELs  in  mind.  Madey  rekindled  the  interest  in  the  FEL  by  his  experiments  in 
the  infrared  wavelength  with  an  rf  linac  at  Stanford  University.24-"25  Around  the 
same  time,  Naval  Research  Laboratory  (NRL)  and  Columbia  University  collab¬ 
orated  on  FEL  experiments  in  the  microwave  regime  with  low  energy  electron 
beams.26-27  In  the  early  1980s,  FEL  experiments  with  high  energy  low  current 
rf  linacs  were  conducted  by  Los  Alamos  National  Laboratory  (LANL),  Spec¬ 
tra  Technology /Boeing  Aerospace  Company  and  TRW/EGG;  FEL  experiments 
with  high  current  low  energy  induction  linacs  were  conducted  at  Lawrence  Liver¬ 
more  National  Laboratory  (LLNL).  Since  that  time  experimental  and  theoretical 
research  rapidly  increased.  Research  on  high  power  FELs  has  been  carried  out  at 
LLNL  and  at  Boeing/LANL.  Past,  current,  and  proposed  experiments  at  LANL, 
TRW,  Spectra  Technology/ Boeing,  Stanford  University,  National  Institute  of 
Standards  and  Technology(NIST)/Naval  Research  Laboratory(NRL),  Duke  Uni¬ 
versity,  Vanderbilt  University,  Rocketdyne/Duke,  Brookhaven  National  Labora¬ 
tory,  Bell  Laboratories,  University  of  California  at  Los  Angeles  (UCLA)  and 
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Lawrence  Berkeley  Laboratory  utilize  higher  energy  beams  to  produce  radiation 
from  the  tenth  of  micron  (in  the  ultraviolet)  to  the  tens  of  micron  (in  the  in¬ 
frared)  wavelength  range.  On  the  other  hand,  University  of  California  at  Santa 
Barbara(UCSB),  Massachusetts  Institute  of  Technology  (MIT),  Columbia,  Uni¬ 
versity  of  Maryland,  University  of  Central  Florida,  Hughes  and  NRL  concentrate 
on  FELs  utilizing  low  energy  beams  to  produce  radiation  from  the  millimeter  to 
the  hundreds  of  microns  in  wavelength.  Currently,  medical  and  material  science 
user  facilities  have  been  established  at  UCSB,  Stanford  University,  NIST/NRL, 
Duke  University,  and  Vanderbilt  University.  The  basic  FEL  theory  has  mostly 
been  verified  with  the  exception  of  optical  guiding.  All  the  experiments  have 
demonstrated  the  importance  of  high  quality  electron  beams.  Research  efforts 
are  underway  to  develop  FELs  in  the  direction  of  user-oriented  and  practical 
devices. 


2  BASIC  CONCEPTS  OF  THE  FEL 


The  basic  principles  of  the  FEL  are  outlined  in  this  section  using  1-D  mod¬ 
els.  We  will  consider  a  linearly  polarized  undulator  with  a  field  which  may  be 
expressed  in  terms  of  the  vector  potential, 

Au(z)  =  Au(z)cos(  f  fcu(z')dz')e*,  (2.1) 

Jo 

where  Au{z)  and  k u(z)  are  the  slowly  varying  amplitude  and  wavenumber  of  the 
undulator.  The  linearly  polarized  radiation  field  is 

Ar(z,t)  =  —A(z)exp[i(kz  —  u>t)\ez/2  +  c.c.,  (2-2) 


where  A(z)  =  |A(z)|e‘^z^  is  the  complex  amplitude  of  the  radiation  field  ex¬ 
pressed  in  polar  variables,  u  is  the  frequency,  k  =  ui/c  is  the  wavenumber,  and 
c.c.  denotes  the  complex  conjugate.  The  electric  field  is  Er  =  —c~1dAr/dt,  and 
the  magnetic  field  is  B  =  Bu  +  Br  =  V  x  (A„  +  Ar). 

The  evolution  of  the  scattered  potentials  is  governed  by  the  wave  equations 


c2  dt2 ) 


Ar 


An  T  . 
-/  —  «/*er, 
c 


(2.3) 


where  Jx  is  the  driving  current  density,  f  —  <Tb/err  is  the  filling  factor  associated 
with  the  radiation  field,  and  &b  and  ar  are  the  cross-sectional  areas  of  the  electron 
beam  and  radiation  beam,  respectively. 

The  transverse  wiggle  and  axial  velocities  are  obtained  by  the  relativistic 
force  equation  for  the  particles, 


dp 

dt 


—  M 


pxJB\ 
7m0c  )  ’ 


(2.4) 
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where  p  =  7 m0v,  7  =  (1  +  p  •  p/mjjC2)1/2  =  (1  —  v  •  v/c2)-1/2  is  the  total 
relativistic  mass  factor,  and  mo  is  the  electron  rest  mass. 

Since  the  Hamiltonian  for  (2.4)  is  not  an  explicit  function  of  transverse 
coordinates,  conservation  of  canonical  transverse  momentum,  i.e., 

7moV.|_  —  je|(Ar  -f  Au)/c  =  constant ,  implies  that 


Vu/C  =  - 


1  |e| 


7  m0c‘ 


(2.5a) 


,  1  |e| 

Vr/C  = - 5 

7  moc2 


(2.56) 


Figure  2  is  a  schematic  drawing  of  an  electron  trajectory  in  the  presence  of  the 
undulator  field.  In  most  cases  of  interest  |AU|  >>  |A|,  in  which  case  we  can 
write 


where  vz  is  the  rms  value  of  the  axial  velocity. 

The  energy  of  an  individual  electron  in  the  electron  beam  changes  according 

~dt  =  ‘ Er)  =  di *’"**>*)•  (2J) 


Since  the  electron  wiggle  velocity  vu  and  the  electric  field  of  the  radiation  Er  are 
in  the  same  direction,  the  electrons  in  the  beam  can  either  gain  or  lose  energy. 
The  electron  dynamics  can  be  simply  described  by  its  motion  in  a  potential 
wave,  called  the  ponderomotive  potential  $ptmrf>28-2B  where 

*  ,  ^  \e\Au-Ar  \e\AuAe^k+k^~^  . 

*pond{z,  t)  = - j—  ^ - 5" - n -  +  C-C-  (2-8 

7  m0cz  7  m0cz  2 

is  formed  by  the  beating  of  the  radiation  field  with  the  undulator  field.  In  writing 
out  the  last  expression  for  the  ponderomotive  potential  the  non-resonant  terms 
have  dropped  from  the  right-hand  side  of  Eq.  (2.8). 

The  oscillating  electrons  at  different  positions  in  the  beam  interact  coher¬ 
ently  with  the  radiation  field  if  their  velocity  vzo  is  approximately  equal  to  the 
phase  velocity 

vptl  =  u>/(k  +  l c„)  (2.9) 

of  the  ponderomotive  wave.  The  resonant  frequency  for  which  vzo  =  is 


u>q  — 


(1  4-  Vzv/c)jQVzokt 

0  +  *2) 


(2.10) 
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where  K  =  \e\Au / y/2mnc2  is  the  normalized  rmd  undulator  vector  potential 
amplitude,  vx0  is  the  initial  rms  axial  velocity,  70  =  7*o(l  +  K2)1^2  is  the  initial 
relativistic  factor,  and  7*0  =  (1  —  v20 /c2)-1^2-  The  FEL  is  tunable  by  varying 
the  magnetic  field  of  the  undulator,  the  period  of  the  undulator,  and/or  the 
energy  of  the  electron  beam. 

Changing  to  Lagrangian  variables  7  and  vz ,  the  equations  of  motion,  as  a 
function  of  z  for  an  electron  with  initial  phase  ij)0  in  the  ponderotive  potential 
well,  are  described  by 


**/c  .  .  .  . 

~T  = - —JbK  a  sin  v  -  eDC 

dz  7 


and 


where 


d?i/> 

~d* 


—  K\  [sin^R  —  sin^] , 


i>  =  f  [(Wc)  +  -  v/vz]  dz'  +  <t>  +  V»o 

Jo 

is  the  phase  of  the  electron  in  the  ponderomotive  wave, 


K  -Ik  (  A 

K-  ~  2K  V(TT^j  j 


1/2 


is  the  synchrotron  wavenumber, 


dk.  1  dK 2  ,  1  +  K 2 

[  dz  ~  2i!  dz  ~  k  7s  eDC 


(2.H) 


(2.12) 

(2.13) 


(2.14) 


(2.15) 


V’r  is  the  resonant  phase,  epc  is  an  applied  DC  accelerating  electric  field,  a  = 
\e\A/ y/2moC2 ,  /r  =  Jo(b)  —  Ji(b)  for  linearly  polarized  undulator,  Jn  is  the  nth- 
order  Bessel  function,  and  b  =  K2 / 2(1  -f  K2).  The  phase  of  the  electron  residing 
at  the  resonant  phase  V'R  does  not  oscillate. 

The  FEL  efficiency,  the  ratio  of  electron  energy  converted  to  radiation  over 
the  initial  electron  energy,  may  be  increased  above  its  intrinsic  value  through 
various  schemes.29-32  Efficiency  can  be  increased  by  increasing  as  a  function 
of  the  axial  distance  z,  the  wavenumber  fcu,  decreasing  the  amplitude  of  the 
vector  potential  K ,  and/or  applying  dc  accelerating  electric  field,  i.e.,  tapering 
such  that  1  >  sinV'R  >  0.  The  total  energy  of  the  electron  is  a  function  of 
the  transverse  wiggle  motion  and  the  mean  axial  motion  such  that  7  =  (1  + 
K2)1^2 7,.  Increasing  the  wavenumber  ku  decreases  the  phase  velocity  of  the 
ponderomotive  wave,  Eq.  (2.9).  The  trapped  electrons  in  the  ponderomotive 
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wave  will  be  decelerated,  thus  kinetic  energy  associated  with  axial  motion  is 
converted  to  radiation.  Decreasing  the  amplitude  of  the  vector  potential  K 
converts  kinetic  energy  associated  with  the  wiggling  motion  to  radiation,  since 
vu,rmt/c  =  Kj 7.  Application  of  a  dc  accelerating  electric  field  eoc  converts 
energy  from  the  dc  electric  field  to  radiation. 

When  the  taper  is  too  fast,  i.e.,  sin  if>R  >  1,  the  electrons  become  free 
streaming  and  no  longer  trapped  in  the  ponderomotive  potential  well.  The 
efficiency  is  further  explained,  and  the  phase-space  diagram  is  described,  in 
Section  3.2. 

The  equation  governing  the  amplitude  of  the  radiation  is  dependent  on  the 
electron  motion33 


da 

dz 


ffBK(z) 


7  /  ’ 


(2.16) 


where  (...)  is  the  ensemble  average  over  all  electrons  initially  within  a  pondero¬ 
motive  potential  wave,  u>&  =  (47r|ej2no/m0)1/2  is  the  plasma  frequency,  and  n0 
is  the  electron  density.  The  derivation  of  (2.16)  from  (2.3)  assumes  the  slowly 
varying  approximation  of  the  radiation  amplitude. 

The  FEL  interaction  in  1-D  is  self-consistently  described  by  Eqs.  (2.11)- 
(2.16).  Three  of  the  following  sections  will  discuss  how  various  nonideal  effects, 
such  as  effective  axial  energy  spread  of  the  electron  beam,  radiation  sideband 
generation,  and  undulator  field  errors,  can  affect  the  performance  of  the  FEL. 
Section  6  shows  how  radiation  guiding  can  improve  the  performance  of  the  FEL. 


3  ELECTRON  BEAM  QUALITY 


One  of  the  most  important  factors  affecting  the  performance  of  the  FEL  is 
the  electron  beam  quality.  The  critical  parameters  determining  the  beam  quality 
for  the  FEL  are  the  peak  current  and  the  effective  energy  spread.  The  shorter 
the  radiation  wavelength,  the  more  stringent  are  the  conditions  on  beam  quality. 
The  efficiency  and  the  gain  (or  growth  rate)  are  measures  of  the  performance  of 
the  FEL.  Table  I  gives  the  expressions  for  gain  and  intrinsic  power  efficiency  in 
different  operational  regimes  for  a  cold  beam.34-35  The  gain  and  efficiency  will 
decrease  as  the  effective  energy  spread  increases. 

If  all  the  electrons  in  the  beam  have  the  same  initial  kinetic  energy,  and 
propagate  along  the  axis  of  the  undulator,  the  electron  beam  can  be  trapped 
in  the  ponderomotive  potential  well  formed  by  an  appropriate  radiation  field. 
The  trapped  electrons  on  the  average  lose  kinetic  energy,  which  is  converted  to 
radiation.  The  largest  efficiency  for  a  given  uniform  undulator  can  be  obtained 
when  the  beam  is  cold. 

For  real  beams,  the  axial  velocities  of  the  electrons  will  have  a  spread.  Poor 
electron  beam  quality  degrades  the  trapping  of  the  electrons  in  the  pondero¬ 
motive  potential  well.  The  FEL  mechanism  fails  when  a  substantial  number  of 
electrons  overtake  and/or  lag  behind  the  ponderomotive  wave. 
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The  source  of  some  of  the  effective  energy  spread  is  related  to  the  undulator. 
An  experimentally  realizable  undulator  field,  which  satisfies  V  x  Bu  =  0  and 
V  •  Bu  =  0,  is  given  by 


Au(x,y,z)  =  Ax(x,y,z)  cos( 


(3.1) 


where  Ax(x,y,z)  takes  on  different  forms  depending  on  the  shape  of  the  mag¬ 
netic  pole  face.  For  pole  faces  with  no  gradient  in  the  sc-direction,  Ax(x,y,z)  = 
Au(z)  cosh(kuy),  where  Au{z)  and  A:u(z)  are  slowly  varying  functions  of  z.  This 
undulator  has  only  weak  focusing  in  the  y-direction.  Here,  we  are  not  consid¬ 
ering  undulators  with  field  errors  (see  Sec.  5)  or  undulators  with  weak  focusing 
in  the  *-direction. 


3.1  Sources  of  Effective  Energy  Spread 


In  the  FEL,  the  electron  beam  can  be  taken  to  be  cold  if  the  rms  axial 
velocity  of  all  the  electrons  in  the  undulator  field  is  identical.  The  major  factors 
that  contribute  to  a  spread  in  vz  are  i)  intrinsic  energy  spread  of  the  beam 
from  the  accelerator,  ii)  transverse  emittance,  iii)  transverse  gradients  of  the 
undulator,  and  iv)  undulator  field  errors.  A  spread  in  vz  is  related  to  a  spread 
in  7 z, 


At* 

7*o 


2  Avz 

7z0 - , 

c 


(3.2) 


where  70  =  (1  +  if2)1/27*o  and  7,  =  7,0  +  At*  =  [l  -  ( vt0  +  A v*)2/c2] 

The  contribution  to  the  total  effective  energy  spread  by  each  of  the  factors  above 
is  discussed  below. 


3.1.1  Intrinsic  energy  spread 

The  intrinsic  energy  spread  of  the  accelerator  can  come  from  the  cathode, 
the  accelerating  structure,  or  the  self-fields  of  the  electron  beam.  The  type  of 
accelerator  affects  the  source  of  the  intrinsic  energy  spread. 

The  electron  beam  can  aquire  an  axial  energy  spread  starting  at  the  cathode.3 
The  causes  for  the  spread  include  i)  temperature  of  the  cathode,  ii)  roughness 
of  the  cathode  surface,  iii)  nonuniformity  of  emission  from  the  cathode  surface, 
iv)  asymmetry,  nonuniformity,  or  nonadiabaticity  of  the  applied  electric  and 
magnetic  fields,  and  v)  the  self-field  of  the  electron  beam. 

The  self-field  of  the  beam  introduces  a  difference  in  energy  across  the  trans¬ 
verse  beam  profile.  The  maximum  energy  spread  of  the  self-field  is  the  difference 
between  the  potential  of  the  beam  on  axis  and  that  on  the  edge  of  the  beam, 

A£,p  =  |eU<Mr6)-<M0)].  (3-3) 
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The  self-potential  is  governed  by  the  equation 


=  — 4ir|e|n0, 


(3.4) 


where  n0  is  the  electron  density.  Substituting  the  solution  of  <£fJ)  from  (3.4)  into 
(3.3),  we  obtain  the  contribution  of  self-fields  to  the  intrinsic  energy  spread 


A  Etp  _  v_ 
E  ~  7o’ 


(3.5) 


where  Budker’s  parameter  is  given  by  v  =  w2r2/4c2  =  /{,/  17/30,  h  is  peak 

current  in  kA,  and  f30  =  (1  —  7 o”2)1^2-  The  contribution  of  self-fields  to  the 
total  energy  spread  is  usually  not  significant  except  for  high  current,  low  energy 
beams. 

For  a  well  designed  cathode,  the  main  source  of  variation  in  the  intrinsic 
energy  spread  among  the  different  accelerators  is  the  accelerating  structure. 
Electrostatic  accelerators,  linear  induction  accelerators,  microtrons,  and  storage 
rings  can  all  have  an  energy  spread  much  smaller  than  1%.  Currently,  rf  linacs 
typically  have  an  energy  spread  of  approximately  1%  due  to  several  possible 
contributions  including  the  varying  accelerating  rf  field  along  the  length  of  the 
electron  bunch. 


3.1.2  Emittance 

Transverse  emittance,38  e,  is  an  important  property  of  beam  quality.  We 
define  7re  to  be  the  area  of  the  smallest  ellipse  in  x,x'  =  dx/dz  transverse  phase 
space  which  encloses  the  particles  in  the  beam.  It  can  be  expressed  approxi¬ 
mately  as 

e  —  rb(0),  (3.6) 

where  is  the  beam  edge  radius  and  6  is  the  ratio  of  the  perpendicular  to  the 
parallel  particle  velocity  in  the  beam.  In  writing  (3.6)  we  have  assumed  that  the 
beam  is  symmetric  and  e  =  tx  =  ey.  In  some  accelerators  the  emittances  in  the 
x  and  in  the  y  planes  are  different.  We  will  consider  the  “envelope,”  or  “edge,” 
emittance  of  the  beam 

eE  =  4  [(*2){(i')J)  -  ((**'))2]1/a  =  4elirm..  (3.7) 

This  definition  is  a  factor  of  four  larger  than  the  rms  emittance. 

As  energy  increases,  e  decreases;  but  the  normalized  emittance 

tn  =  /?7  e  (3.8) 

is  invariant  through  the  entire  linac  and  beam  transport  system  for  linear  fo¬ 
cusing  systems.  The  FEL  requirement  on  the  smallness  of  en  is  therefore  a 
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specification  of  the  injector  system.  With  new  reseach  on  laser-initiated  photo- 
emission  cathodes,  future  linacs  have  the  potential  of  reaching  higher  currents 
at  low  emittance. 

A  frequently  used  figure  of  merit  for  beam  quality  in  FELs  is  the  normalized 
beam  brightness 

Bn  =  a0Ib/n2e7n ,  (3.9) 

where  Iy  is  the  peak  beam  current  and  ao  is  a  dimensionless  number  of  order 
unity.  For  beams  with  a  uniform  ellipsoidal  cross  section,  ao  =  2.  Table  II  lists 
beam  parameters  associated  with  various  accelerators  used  for  FELs. 

Finite  transverse  emitldnce  affects  FELs  in  two  ways:  i)  it  introduces  an 
additional  source  of  effective  energy  spread,  and  ii)  it  causes  divergence  of  the 
electron  beam  which  limits  the  radiation  wavelength. 

Consider  an  electron  beam  propagating  in  the  linearly  polarized  undulator 
defined  by  Eq.  (3.1).  The  undulator  has  a  transverse  gradient  in  the  y-direction 
only.  Since  there  is  no  focusing  mechanism  in  the  x-direction,  the  beam  diverges. 
If  the  beam  is  originally  cold,  i.e.,  all  the  electrons  have  identical  energy,  the 
beam  divergence  in  the  x-direction  causes  an  effective  spread  of  the  velocity  in 
the  z-direction,  vZith, 


_ TzO _ 


-7.0  |l-7.!o(*i  +  ^)/2]- 


(3.10) 

Consequently,  the  effective  energy  spread  due  to  emittance  in  the  x-direction  is 


7.o  “2(1  +  K*)\rJ  • 


(3.11) 


The  contribution  of  emittance  to  the  total  effective  axial  energy  spread  is  often 
the  limiting  factor  in  applying  an  accelerator  to  the  generation  of  short  wave¬ 
length  radiation  with  the  FEL  mechanism. 

Furthermore,  emittance  considerations  lead  to  an  upper  bound  on  the  ra¬ 
diation  wavelength  based  on  geometrical  overlap  of  the  electron  beam  and  the 
radiation  beam.3  In  the  x-direction,  the  electron  beam  profile  diverges  as 


rk(z)  =  rt(0)|l+eV/r<(0)],/!,  (3.12) 

where  rb( 0)  is  the  electron-beam  radius  at  z  =  0.  For  a  Gaussian  radiation 
mode,  the  radial  extent  of  the  beam  is 

r.(z)  =  M0)[l+A2z2/w2rJ(0)]1/2,  (3.13) 
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where  r,(0)  is  the  radiation  -beam  waist  and  A  is  the  radiation  wavelength.  If  we 
require  rj,  <  r,  to  ensure  good  geometrical  overlap  of  the  radiation  and  electron 
beams,  we  find  the  emittance  requirement  to  be 

A  >,«.  (3.14) 

For  electron  beams  with  weak  focusing  in  both  x-  and  y-directions,  this  in¬ 
equality  is  not  applicable;  also  this  inequality  is  not  a  sufficient  condition  for  the 
beam  quality  requirement  in  an  FEL  because  of  the  contribution  of  emittance 
to  the  effective  energy  spread. 

3.1.3  Weak  focusing  of  the  undulator 


In  the  y-direction,  the  weak  focusing  of  the  transverse  gradient  (3.1)  of 
the  undulator  field  causes  the  electrons  to  undergo  betatron  oscillations.  The 
betatron  oscillations  introduce  an  effective  energy  spread37-39  even  for  cold 
beams  with  zero  emittance.  The  derivation  is  given  below. 

The  particle  motion  in  the  y-direction  is  governed  by 


dPy 

dt 


(3.15) 


where  Bz  =  —  Auku  sinh(fcuy)  cos(fcuz)  is  the  z-component  of  the  magnetic 
undulator  field  associated  with  the  realistic  undulator  and  t5*  =  vxo  + 

V2K  cosh(/cuy)  cos(fcuz).  Here,  the  particle  orbits  y  and  vx  are  functions  of 
(z,  Xo ,  yo ,  vx0i  VyQ ),  where  d/dt  =  vz(d/dz),  x0  and  y0  are  the  initial  transverse 
coordinates,  and  vx0  and  vy0  are  the  initial  transverse  velocities.  Assuming  the 
fast  oscillatory  terms,  with  wavelength  equal  to  half  the  undulator  wavelength, 
as  unimportant  and  replacing  vz  by  c  at  appropriate  places,  we  find  that 

+  !r~2ku  sinh(2fcuy)  =  0.  (3.16) 


The  small  amplitude  solution  for  motion  in  the  y-direction  is 

V  —  yb  cos(Kgz  -  d>p), 


(3.17) 


where  Kg  —  Kku/ 70  is  the  wavenumber  of  the  betatron  oscillations, 


is  the  amplitude  of  the  betatron  oscillation,  and 

4>g  =  cos-I(y0/y(,) 


(3.18) 


(3.19) 
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is  the  initial  phase  of  the  betatron  oscillation. 

Applying  conservation  of  energy  in  a  static  magnetic  undulator  field,  we 
can  show  that  the  axial  velocity  of  the  particles  decreases  as  the  amplitude  of 
the  betatron  oscillations  increases.  This  results  in  an  equivalent  energy  spread. 
The  axial  velocity  is  found  to  be 

vl  7o  -  1 


- 2~  cosh2  (kuy)  cos2  (kuz)  +  ~[cosh(2fc1ty0)  —  cosh(2feuy)] - y. 

To  *  _  £ 

(3.20) 

Dropping  terms  that  oscillate  at  twice  the  undulator  wavenumber,  we  obtain 
vz  =  vz0  —  Auz>0,  where  vto  =  foe  —  (K2/2jg)c  is  the  mean  axial  velocity  of  the 
electron  travelling  on  axis,  /30  =  (1  -  T^-2)1^2,  and 


2 


(3.21) 


is  the  amount  that  the  axial  velocity  is  reduced  for  particles  executing  betatron 
oscillations  with  oscillation  amplitude  r&  =  y&.  When  the  electron  beam  is 
properly  matched  into  the  undulation,  i.e.,  the  beam  envelope  does  not  oscillate, 
then  the  beam  radius  is  related  to  the  emittance 


rl  =  c/vKp. 


3.1.4  Total  effective  axial  energy  spread 
An  estimate  of  the  total  effective  axial  energy  spread  is 


(3.22) 

The  terms  on  the  right-hand  side  of  Eq.  (3.22)  represent  the  contributions  of 
the  effective  energy  spread  from  the  intrinsic  energy  spread  of  the  accelerator, 
transverse  emittance  in  the  x-direction,  transverse  emittance  associated  with 
weak  focusing  of  the  undulator  in  the  y-direction,  and  the  undulator  field  error. 
The  discussion  relating  to  the  undulator  field  errors  will  be  presented  in  Section 
6.  The  estimate  of  the  total  effective  axial  energy  spread  (3.22)  is  approximate 
because  we  have  assumed  that  the  various  contributions  to  it  are  statistically 
independent. 


3.2  Acceptable  Energy  Spread  Limit 

Estimates  of  acceptable  limits  of  energy  spread  can  be  obtained  using  sim¬ 
ple  trapping  arguments.28-29’84*35  In  the  linear  development  of  the  laser  radia¬ 
tion  the  injected  axial  beam  velocity  is  slightly  greater  than  the  phase  velocity, 
w*o  =  Vph  +  Aw*,  where  vph  =  u>/(fc  +  fcu)  is  the  axial  phase  velocity  of  the  longi¬ 
tudinal  wave,  vph  >>  Aw*  >  0,  and  Avz  depends  on  the  particular  FEL  regime 
under  consideration.  This  is  illustrated  in  Fig.  3a,  where  fd(~l)  is  the  electron  dis¬ 
tribution  function  at  the  entrance  to  the  undulator  and  7pk  =  (1  —  v2h/c2)-1/2. 
The  radiation  amplitude  increases  at  the  expense  of  the  electron’s  kinetic  en¬ 
ergy  until  the  electrons  become  deeply  trapped  in  the  longitudinal  wave.  At  this 
point  the  radiation  field  reaches  its  maximum  amplitude,  and  the  average  axial 
electron  velocity  is  approximately  given  by 

Vz\sat  =  Vph  -  Avz.  (3.23) 

The  electron  distribution  function  at  the  exit  of  th;  idulator  is  illustrated 
in  Fig.  3b.  At  saturation  the  average  axial  electron  velocity  has  decreased  by 
approximately  2 Aw*.  The  decrease  in  the  electron  beam  energy  can  be  directly 
equated  to  the  increase  in  radiation  energy.  For  highly  relativistic  electron  beams 
the  decrease  in  the  average  electron  kinetic  energy  is  A E  =  2'y0'yl0movz0Avz 
and  hence  the  radiation  efficiency  is 

A  F 

"  =  5,— h;?  =  (324) 

The  longitudinal  wave  “sees”  the  beam  as  monoenergetic  if  the  beam’s  axial 
velocity  spread  is  small  compared  to  Aw*.  Requiring  Aw*>tk  <<  2Aw*,  we  obtain 


TzO 


<  V- 


(3.25) 


One  way  to  visualize  the  trapping  condition  is  to  plot  the  electron  distri¬ 
bution  in  the  phase  space  ( xl> ,  dip / dz ),  where  ^  is  the  relative  phase  between  the 
electrons  and  the  ponderomotive  wave.  It  is  more  convenient  to  visualize  the 
phase  space  (Vst)  or  (Vs7z)>  since  the  phase  xp  is  related  to  the  parallel  energy 
7*  by 

diji/dz  =  2fcu(7*  -  7ro)/7zo-  (3.26) 

Figure  4  is  a  schematic  of  the  phase  space  (Vs7z)  for  uniform  undulators.  The 
separatrix,  which  separates  trapped  and  flowing  orbits  in  the  ponderomotive 
potential,  is  outlined.  The  bucket  width  changes  adiabatically  as  the  radiation 
amplitude  grows.  The  electrons  become  trapped  inside  the  separatrix.  When  the 
electrons  rotate  to  the  lowest  point  of  the  separatrix,  the  FEL  attains  the  largest 
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intrinsic  efficiency.  The  bucket  width  is  approximately  the  intrinsic  efficiency. 
This  is  a  graphic  illustration  of  the  inequality  requirement  in  Eq.  (3.25). 

In  physical  terms,  the  effective  energy  spread  of  the  beam  does  not  degrade 
the  performance  of  the  EEL  if  electrons  with  different  axial  velocities  do  not 
significantly  separate  at  the  end  of  some  relevant  interaction  length.  This  sepa¬ 
ration  distance  must  be  less  than  some  fraction  (dependent  on  the  FEL  regime) 
of  the  ponderomotive  potential  wavelength.  The  wavelength  of  the  ponderomo- 
tive  potential  is  2ir/(k  +  ku),  which  is  approximately  the  radiation  wavelength 
for  large  7z0.  The  effective  energy  spread  criterion  can  be  written  as 

(fcu  +  k)AvZtthLtff/vz0  «  2 ttF,  (3.27) 

where  A vZith  is  the  axial  velocity  spread,  ie//  is  the  effective  interaction  length, 
and  F  is  the  fraction  of  the  2ir  phase  the  electrons  would  have  traveled  in  the 
ponderomotive  potential  in  the  effective  distance  for  a  cold  beam. 

In  the  low  gain  regime, 

Lz //  =  Lu,  (3.28) 

and  the  the  electrons  obtain  the  highest  efficiency  when  they  traverse  almost  the 
whole  2ir  phase  before  bunching  at  the  bottom  of  the  pondermomotive  potential 
well,  i.e.,  F  ~  1.  Substituting  Lu  into  Eq.  (3.27),  we  obtain 


A  7*,tft. 
7*o 


=  Vi 


(3.29) 


where  N  is  the  number  of  undulator  periods  and  rj  is  the  intrinsic  efficiency. 

In  the  high  gain  Compton  regime,  the  electron  density  is  sufficiently  low  so 
that  collective  space-charge  effects  are  neglegible.  In  this  regime,40 


Leff  —  Le, 


(3.30) 


where  Le  =  l/T  is  the  e-folding  length,  and  T  is  the  peak  growth  rate.  For  the 
beam  to  behave  as  cold,  the  electrons  should  not  deviate  by  more  than  n/2  in 
the  27r  phase  of  the  ponderomotive  well,  i.e.,  F  ~  1/4.  Substituting  Le  into  Eq. 
(3.27),  we  obtain 

(3.31) 

7*o 

In  the  Raman  regime,  the  force  on  the  beam  electrons  due  to  the  collec¬ 
tive  space-charge  wave  becomes  comparable  to  the  ponderomotive  wave.  We 
can  estimate  the  effective  interaction  length  to  be  the  Doppler-shifted  plasma 
wavelength 

Leff  -  vz07zOi/yo/vb,  (3.32) 
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where  is  the  nonrelativistic  plasma  frequency.  Taking  F  ~  1/4,  we  obtain 


^^<<77.  (3.33) 

7*o 

If  the  intrinsic  efficiency  is  small,  one  can  tolerate  a  larger  thermal  spread 
in  the  beam  by  imposing  a  large  radiation  field  to  initially  trap  the  electrons  and 
extract  their  kinetic  energy  by  appropriately  tapering  the  undulator.  A  large 
fraction  of  the  electrons  in  the  ponderomotive  potential  well  can  be  trapped,  and 
their  energy  can  be  extracted  by  adiabatically  tapering  the  undulator.  Three 
methods  (also  known  as  efficiency  enhancement  methods)  can  be  applied:  de¬ 
crease  the  period  of  the  undulator,  decrease  the  vector  potential  of  the  undulator, 
and/or  apply  an  accelerating  electric  field. 

Figure  5  plots  the  electron  distribution  function  at  the  entrance  and  the  exit 
of  the  undulator.  To  trap  the  electrons  initially,  we  require  (70)  =  7h(0),  where 
7R(*)  =  [l+*-M’]‘/J7,kW  is  the  7  associated  with  the  resonant  electron  (see 
Fig.  5a).  At  the  end  of  the  undulator,  the  trapped  electrons  have  lost  energy 
(see  Fig.  5b). 

Figure  6  is  a  plot  of  the  electrons  in  the  phase  space  (^,7)  and  the  separatrix 
for  the  resonant  phase  ipR  that  satisfies  0  <  sin^H  <  1.  To  trap  a  substantial 
fraction  of  the  electrons,  we  require 

7fl(0)  -  A-ytrap/l  <  7  <  7*(0)  +  A-ytraP /2,  (3.34) 

where  A^trap  —  A7Iitrap[l  -f  Ar2(z)]1'/2  is  the  full-width  of  the  bucket,  as  derived 
below.  Electrons  initially  trapped  within  the  separatrix  remain  trapped  provided 
they  are  adiabatically  decelerated.  The  half-width  of  the  bucket  is  given  by 

=±2  K.giM,  (3-35) 

where  K ,  is  the  synchrotron  wavenumber  (2.14)  and 

s(V’r)  =  [(V>R  -  *-/2)sin^R  +  cosV’r]172  (3.36) 

is  the  reduction  in  the  width  of  the  bucket  as  a  function  of  the  rate  of  taper.  In 
order  to  trap  a  substantial  fraction  of  the  electrons,  we  require  the  inital  depth  of 
the  trapping  potential  |e|<?!>ZitT.0p  =  A7 z<ira.p{'ii<s / 7 zo)m0c7  to  be  larger  than  the  ef¬ 
fective  axial  energy  spread,  i.e.,  \e\(f>z  ,trap  ^  5”  AEth,total  —  A~yz<  1/1(70  / H zO )rn<>c  . 
By  combining  (3.26)  and  (3.35),  the  full-width  of  the  trapping  bucket  is 

A7z,trap  _  _1_  dip 
7*0  fcu  dz 


=  4  g{ii>R) 


(fsK  lalV 

\1  +  K*J 


(3.37) 
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According  to  Eq.  (3.37),  even  when  condition  (3.34)  is  satisfied,  the  fraction  of 
trapped  electrons  decreases  as  the  taper  increases  because  the  phase-space  area 
decreases.31 

The  effective  axial  energy  spread  of  the  electron  beam  is  unavoidable  be¬ 
cause  it  can  be  introduced  through  the  self-fields,  transverse  finite  emittance, 
and  transverse  gradients  of  the  undulator.  The  above  expressions  specify  the  cri¬ 
teria  that  this  effective  axial  energy  spread  must  satisfy  for  uniform  undulators 
and  for  tapered  undulators. 

Electron  beams  can  possess  other  energy  variations  that  are  also  undesir¬ 
able.  The  beam  can  have  an  energy  jitter  or  an  energy  droop  as  a  function 
of  time  arising  from  the  specific  accelerator  design.  If  the  electron  beam  is 
driven  by  an  rf  source,  the  micro-pulses  can  have  phase  jitter  relative  to  the  rf 
source.  The  beam  energy  may  vary  from  one  macro-pulse  to  another.  All  these 
non-ideal  effects  degrade  the  performance  of  the  FEL. 

4  RADIATION  SIDEBANDS 

A  major  problem  affecting  FEL  performance  is  the  growth  of  sideband  radi¬ 
ation.  These  sidebands  are  shifted  in  frequency  from  the  fundamental  radiation 
and,  as  they  grow,  they  deplete  the  electron  beam  of  energy,  thus  reducing 
the  growth  of  the  fundamental.  Sidebands  in  FELs  arise  as  follows.  Electrons 
trapped  in  the  ponderomotive  well  perform  synchrotron  oscillations.  Particles 
trapped  near  the  bottom  of  the  well  have  a  simple  harmonic  motion  with  period 
given  by  2iv/Kt,  where  K,  is  the  synchrotron  wavenumber,  Eq.  (2.14).  The 
oscillations  of  the  trapped  electrons  generate  current  components  at  frequencies 
shifted  from  the  electron  wiggling  frequency  cku  by  some  multiple  of  the  syn¬ 
chrotron  frequency  cK,.  The  synchrotron  radiation  will  appear  at  frequencies 
separated  from  the  desired  carrier  frequency  by  a  multiple  of  twice  the  Doppler 
upshifted  synchrotron  frequency,  i.e., 

A  u>,  =  cA  k,,  (4.1) 


where  A k,  =  2~flK./{l  +  K2)  and  K,  =  2 ku  [fBK\a\/(l  +  A2)]1/2.  This  growth 
of  radiation  sidebands,  referred  to  as  the  sideband  instability,  was  first  pointed 
out  by  Kroll  et  al.31 

The  analysis  by  Kroll  et  al.  revealed  the  physical  basis  for  this  instability  to 
be  similar  to  that  in  stimulated  Raman  scattering.  The  signal  wave  u>  can  scatter 
into  an  anti-Stokes  (higher  frequency)  wave  u;'  by  extracting  energy  from  the 
synchrotron  motion  of  the  electrons  (u>  4-  cAfc,  — ♦  u>’).  Alternatively  the  signal 
wave  can  scatter  into  a  Stokes  (lower  frequency)  component  u>'  andt  at  the 
same  time  cause  the  electrons  to  go  into  higher  energy  states  of  the  synchrotron 
motion  (u;  — ♦  uj'  +  cAfc,).  The  growth  rate  and  the  relative  amplitude  of  these 
two  processes  is  determined  by  the  distribution  function  of  the  electrons  in  the 
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vicinity  of  vph.  The  Stokes  process  is  obviously  the  more  worrisome  of  the  two 
since  it  can  lead  to  detrapping  of  the  electrons  out  the  ponderomotive  buckets. 

Since  the  work  of  Kroll  et  al.  more  detailed  theoretical  analyses,41-45  nu¬ 
merical  simulations, 48-57  and  experimental  observations45,53,58  of  the  sideband 
instabilities  have  been  reported  in  the  literature. 

The  sideband  instability  may  be  important  when  the  electrons  can  make 
approximately  one  or  more  synchrotron  oscillations  in  the  undulator.  This  is 
the  case  for  an  FEL  amplifier  after  saturation,  an  FEL  operating  in  the  trapped 
particle  mode,  or  an  FEL  oscillator  with  high  intra-cavity  power. 

The  synchrotron  frequency  is  a  function  of  the  signed  wave  amplitude  |a|. 
Hence,  in  a  high  gain  FEL  where  |a|  changes  rapidly  in  a  synchrotron  wavelength, 
the  resonance  condition  is  met  only  for  a  short  time  and  the  instability  is  unlikely 
to  be  important.  For  an  oscillator,  the  radiation  intensity  inside  the  resonator 
is  substantially  higher  than  that  of  the  extracted  radiation,  and  the  amplitude 
slowly  evolves  between  the  cavity  mirrors.  This  situation  tends  to  produce 
significant  levels  of  sidebands  and  may  possibly  lead  to  diffusion,  or  Brownian 
motion,  of  the  electrons  out  of  the  ponderomotive  well.31 

Detrapping  is  obviously  a  serious  effect  since  it  leads  to  a  reduction  in 
gain.  Additionally,  the  sideband  instability  modulates  the  ouput  signal  and 
consequently  increases  its  spectral  width.  The  performance  of  the  mirrors  in 
an  oscillator-type  operation  can  be  harmed  from  the  modulation  of  the  wave 
envelope  caused  by  the  sidebands. 

A  schematic  illustration  of  radiation  with  a  large  sideband  component  is 
shown  in  Fig.  7.  The  amplitude  of  the  electric  field  is  modulated,  as  shown  in 
Fig.  7a.  The  wavelength  of  the  modulation  is  approximately  A,  =  2n/Akt.  The 
oscillating  radiation  amplitude  leads  to  an  oscillating  separatrix,  as  shown  in 
Fig.  7b.  The  electrons  near  the  border  of  the  separatrix  gradually  become  lost 
as  the  separatrix  oscillates.  The  loss  of  trapped  electrons  results  in  the  loss  of 
efficiency. 

For  an  FEL  amplifier,  sideband  modulation  can  be  reduced  by  a  rapid 
tapering  of  the  undulator  field  parameters.45,55’57  In  a  rapidly  tapered  FEL 
amplifier,  the  ponderomotive  bucket  changes  its  position  in  phase  space  fast 
enough  to  distort  the  synchrotron  orbits  of  the  electrons.  A  measure  of  the 
distortion  of  electron  orbits  is  given  by31,55’57 


cdm/dz 
(cK.{ 7  -  7r)) 


(4.2) 


which  i6  the  ratio  of  the  change  in  energy  of  the  resonant  particle  cd-)R/dz  due  to 
tapering  and  the  change  in  energy  (cK,( 7-7/?))  due  to  synchrotron  motion.  For 
a  tapering  rate  corresponding  to  a  resonant  phase  of  V’fi  =  18°,  the  distortion  of 
electron  orbits  amounts  to  R  =  25%.  Together  with  a  decrease  in  the  trapping 
fraction  at  fast  tapering  rates,52  substantial  reduction  in  sideband  modulation 
was  observed  in  numerical  simulations  of  an  FEL  amplifier  configuration.55,57 
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Several  ideas  have  been  put  forward  in  order  to  reduce  the  sideband  ampli¬ 
tude  in  FELs;  these  include 

1.  Use  of  mirrors  with  poor  reflectivity  at  the  undesirable  sideband  wave- 
lengths.47’51’83-84 

2.  Introduction  of  gratings86  inside  the  oscillator  cavity  so  that  the  sideband 
radiation  is  deflected  outside  the  path  of  the  oscillator. 

3.  Introduction  of  a  gas  or  liquid  with  sufficient  dispersion  to  suppress  the 
build-up  of  the  sideband  in  the  cavity.31 

4.  Rapid  tapering  of  the  undulator.48’55’87 

5.  Operation  in  the  regime  where  diffraction  is  important.89-80  Theoretical 
calculations  indicate  that  in  this  regime  the  sideband  growth  rate  can  be 
reduced  under  certain  conditions,  but  not  eliminated. 

Thus  far  all  analyses  of  the  sideband  instability  are  based  on  highly  sim¬ 
plified  models  of  the  FEL  wherein  many  effects  are  neglected:  the  evolution 
of  the  optical  field  amplitude,  the  transverse  variation  of  the  optical  field,  the 
diffraction  of  the  optical  field,  etc.  Inclusion  of  these  effects  and  a  self-consistent 
determination  of  the  electron  distribution  is  required  before  a  quantitative  com¬ 
parison  between  theory  and  experiment  can  be  attempted. 

5  UNDULATOR  FIELD  ERRORS 

Intrinsic  magnetic  field  errors  SB  are  present  in  any  realistic  undulator  mag¬ 
net.  Such  errors  are  unavoidable  and  arise  from  imperfections  in  the  fabrication 
and  assembly  of  undulator  magnets.  State-of-the-art  undulator  construction 
techniques  yield  rms  field  errors  on  the  order61  (SB /Bu)rm,  —  0.1  —0.5%.  These 
field  errors  perturb  the  electron  beam  as  it  propagates  through  the  undulator 
and  lead  to  i)  a  random  walk  of  the  beam  centroid,62-69  Sx ,  ii)  variations  in 
the  parallel  beam  energy,85-89  Sjz,  and  iii)  variations  in  the  relative  phase  of 
the  electrons  in  the  ponderomotive  potential,85-89  Sip.  If  left  uncorrected,  field 
errors  ultimately  decrease  free-electron  laser  (FEL)  gain62-69  (this  reduction  be¬ 
comes  more  significant  for  long  undulators).  Reduction  in  gain  may  occur  from 
a  loss  of  optical  guiding  (due  to  large  Sx)  or  from  a  loss  of  FEL  resonance  (due 
to  large  Sip). 

Given  the  precise  functional  dependence  of  the  undulator  errors  SB(z)  for 
a  given  undulator,  one  may  obtain  Sx(z)  for  that  specific  undulator.  However, 
one  does  not  always  know  ahead  of  time  the  full  functional  dependence  of  SB(z). 
Instead,  one  may  know  only  certain  statistical  properties  of  the  field  errors,  such 
as  the  rms  value  SBrmt.  Hence,  it  is  useful  to  consider  an  ensemble  of  statisti¬ 
cally  identical  undulators  for  which  the  statistical  properties  of  the  field  errors 
are  known.  By  performing  appropriate  averages  over  this  ensemble,  one  may 
determine  the  mean  (Q)  and  variance  <r  for  a  quantity  Q  and,  hence,  determine 
the  most  probable  range  of  a  single  realization  of  Q.  Throughout  the  following, 
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SB(z)  is  assumed99'89  to  be  a  random  function  with  zero  mean,  and  finite  vari¬ 
ance  and  with  an  autocorrelation  distance  given  by  zac  ~  \u/2.  Also,  in  the 
following,  a  helical  undulator  will  be  assumed,  and  generalization  of  the  results 
for  a  linear  undulator  is  straightforward.99’99 

5.1  Random  Walk  of  the  Beam  Centroid 

As  the  electron  beam  propagates  through  the  undulator,  the  electrons  ex¬ 
perience  random  velocity  kicks  Svx  via  the  vz  x  6 By  random  force.  The  mean- 
square  centroid  motion  for  the  electron  beam,  including  the  effects  of  undulator 
transverse  gradients  (weak  focusing),  is  given  by98,tN 

M  -  !  ("  T?2)  ■ ■  M 

where  D  =  alkl(SBl)zac/(2j2),  kg  =  auku/y/ 2y0  is  the  betatron  wave  number 

of  the  helical  undulator,  SB  =  SB/BU,  au  =  (|e| /moC2)Bu/ku  is  the  normalized 
vector  potential  of  the  helical,  undulator  and  Bu  is  the  ideal  undulator  peak 
magnetic  field.  Physically,  the  centroid  orbits  6x  and  6(3X  =  6vz/c  represent 
diffusing  betatron  orbits  characterized  by  a  diffusion  coefficient  D.  Notice  that 
by  increasing  k by  additional  external  focusing,  one  may,  in  principle,  keep 
Sxrmt  as  small  as  desired.  Furthermore,  notice  that  in  the  1~D  limit,  (2 kgz)2  « 
1,  (Sfll)  =  2Dz,  and  ( Sx 2)  —  2Dz3 / 3.  Hence,  weak  focusing  (finite  kg)  is 
effective  in  reducing  the  asymptotic  scaling  of  the  random  walk  Sxrmt  from  z3/2 
to  z1/2.  To  avoid  loss  of  optical  guiding  it  is  desirable  to  keep  (Sx2)  <<  r], 
where  r,  is  the  radiation  spot  size. 

5.2  Variations  in  the  Parallel  Beam  Energy 

Not  only  do  the  field  errors  perturb  the  perpendicular  motion  of  the  elec¬ 
trons,  they  also  perturb  the  parallel  motion.  This  is  true  since  a  static  mag¬ 
netic  field  conserves  total  electron  energy.  The  parallel  motion  may  easily 
be  calculated88’89  using  the  above  expressions  for  the  perpendicular  motion 
along  with  4-  (32x  =constant.  One  may  calculate  various  statistical  mo¬ 
ments  of  the  parallel  motion,88’99  such  as  the  mean  parallel  energy  variation 

(ht)  =  (lx)  -  7x0, 


(hr) 

7x0 


(l+«S/4) 

(1+a2)’ 


a\k2u(SB2)za 


(5.2) 


where  the  limit  (2kgz)2  >>  1  has  been  assumed  along  with  (SB\)  =  (SB 


2). 

V' 
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Statistically,  (£7*)  may  be  interpreted  as  an  effective  energy  spread  due  to 
field  errors.88’69  This  effective  energy  spread  may  lead  to  loss  of  FEL  resonance. 

5.3  Deviations  in  the  Relative  Phase 

To  determine  how  the  parallel  energy  variation  affects  FEL  gain,  it  is  nec¬ 
essary  to  consider  the  relative  phase  t/*  of  the  electrons  in  the  ponderomotive 
wave,  dij)/dz  =  k  +  Jb„  -  u >/{c0z).  In  the  small  signal  limit  (a  — »  0,  where  a  is 
the  normalized  radiation  field),  the  deviation  in  phase  6il>  due  to  the  field  errors 
is  given  by 

dz'(20xoS/3±  +  k26x2±  +  601),  (5.3) 

where  0XO  is  the  ideal  wiggle  motion  in  the  absence  of  field  errors,  i.e.,  (3Z 0  = 
(au/yo)cos(kvz),  and 

=  fZ  dz'sinksiz' -z)6By(z')  (5.4) 

7*/3  Jo 

and  0T  =  dhx/dz.  Similar  expressions  hold  for  6y  and  6(3y.  Statistically  averag¬ 
ing  over  the  undulator  ensemble  gives 

m--iA (5.5) 

(!  +0 

where  {SB2)  =  {SB2)  =  {SB2)  has  been  assumed.  Notice  that  the  mean  phase 

deviation  (6ij))  is  independent  of  the  effects  of  transverse  focusing. 

Physically,  Sxf}  may  be  interpreted  as  an  oscillation  of  the  ponderomotive 
well  due  to  field  errors.  Maintaining  FEL  resonance  requires  Si}>  to  be  small 
compared  to  tt,  i.e.,  the  width  of  the  well.  In  the  low  gain  regime,  this  phase 
deviation  must  be  kept  small  over  the  entire  undulator  length  Lu.  Requiring 
|  {8^){z  =  Lu ))  |  <  2w  implies  6Brm,  <  a/(irN),  where  a2  =  (1  +  a„)/(2a* ).  A 
similar  condition  may  be  obtained  from  considerations  on  the  effective  energy 
spread,  (£7*)/7xo  <V  =  l/(2N).  In  the  high  gain  regime,  the  situation  is  some¬ 
what  different,  since  the  length  scale  over  which  the  FEL  resonant  interaction 
occurs  is  the  e-folding  length  l/T,  where  T  is  the  spatial  growth  rate  of  the 
radiation.  Maintaining  resonance  in  the  high  gain  regime  corresponds  to  keep¬ 
ing  Sip  small  over  an  e-folding  length:  |  {Sxp{z  =  1/T))  |  <  rr/2.  Since,  typically 
1/r  <<  Lu,  one  expects  the  high  gain  not  to  be  strongly  affected  by  the  phase 
deviation  Sxf)  (in  contrast  to  the  low  gain). 
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5.4  Degradation  of  FEL  Gain 

Quantitatively,  the  effect  of  the  phase  deviation  on  the  FEL  gain  in  the  low 
gain  regime  may  be  determined  analytically.  The  normalized  mean  amplitude 
gain  is  related  to  Sip  by  the  following  expression: 

( G )  =  f  dz'  /  dz"(z'  —  z")  (sin  [nku(z'  —  z")  +  A5^/>]) ,  (5.6) 

Jo  Jo 

where  A  Sip  =  Sip(z')  —  Sip(z")  and  n  =  —NX(uj  —  u>0)/2nc  is  the  normalized 
frequency  mismatch  parameter.  Setting  A Sip  =  0  in  the  above  equation  gives 
the  gain  in  the  absence  of  field  errors. 

Evaluation  of  the  ensemble  average  in  the  above  expression  is  dependent  on 
the  statistical  distribution  of  the  function  A Sip.  Recall  that  the  phase  deviation 
6 ip  is  proportional  to  terms  that  are  linear  in  the  field  error  SB  as  well  as  terms 
that  are  quadratic  in  the  field  error,  as  indicated  by  Eq.  (5.3).  If  the  field 
error  SB  is  a  Gaussian  distributed  random  variable,  then  terms  quadratic  in 
SB  obey  a  gamma  distribution.  Hence,  if  the  quadratic  terms  dominate  in  the 
expression  for  Sip,  then  Sip  will  tend  to  be  gamma  distributed.  Assuming  A  Sip 
to  be  approximately  gamma  distributed  allows  the  ensemble  average  in  Eq.  (5.6) 
to  be  evaluated  using  the  Rice-Mandel  approximation,82’89’70  yielding 

(G)  =  J’  iz'  jf  dz"(z'  -  z")  (1  +  (A (5  7) 

x  sin  [nku(z'  ~  z")  +  f  tan-1  ((A Sip)/ f)]  , 
where  /  =  (A Sip)2 /  ((A Sip2)  —  (A Sip)2). 

It  is  possible  to  show  that  the  mean  gain  is  a  function  of  only  two  param¬ 
eters,  ( G )  =  F  [n,  (^V’)mor),  where  {S^)max  =  (6V>(2  =  i«))*  Furthermore,  one 
can  show  that  the  maximum  mean  gain  ( G)max  decreases  as  increases. 

Equation  (5.7)  may  be  evaluated  numerically  to  determine  the  behavior  of 

the  mean  gain.  Figure  8  illustrates  this  behavior,  in  which  the  mean  gain  ( G ) 
is  plotted  as  a  function  of  the  frequency  mismatch  n  for  several  values  of  rms 
field  error  SBrmt.  The  parameters  in  Fig.  8  correspond  to  a  linearly  polarized 
undulator  with  Bu  =  5.4  kG,  =  2.8  cm,  Lu  =  364  cm,  and  7  =  350  in  the 
limit  k/3  =  0  (transverse  focusing  is  neglected). 

5.5  Beam  Steering 

One  method  for  reducing  the  detrimental  effect  of  field  errors  is  through 
the  use  of  beam  steering.62-89  The  effect  of  steering  is  to  reduce  the  mean  phase 
deviation  by  a  factor  of  1/3,  (Si>(Lu)\  rteering  )  —  ( 1/3)(^V,(  Lu  )  |no  steering)  i  when 
one  steering  section  is  used.  It  is  also  possible  to  calculate  ( G )  including  the 
effects  of  steering. 
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Steering  has  reduced  the  variance  of  the  phase  deviation  by  a  significant 
amount.  For  cases  in  which  kp  ^  0,  it  is  possible  to  show88’69  that  steering 
reduces  the  mean  phase  deviation  when  the  length  over  which  the  steering  is 
performed  is  less  than  the  betatron  wavelength,  Lt  <  \p.  For  cases  in  which 
L.  >  \p,  beam  steering  may  increase  the  value  of  (Sift). 

This  concept  may  be  generalized  to  the  case  of  multiple  beam  steering,  in 
which  the  electron  beam  is  steered  back  to  the  axis  in  several  places  along  the 
length  of  the  undulator.83-84’68-89 

Self-consistent,  three-dimensional,  numerical  simulations  with  undulator 
field  errors  and  beam  steering  were  performed  for  a  linearly  polarized  undulator. 
The  amplitude  of  the  undulator  field  is  assumed  to  be  in  error  with  the  error 
of  the  form  of  half  a  period  of  a  sinusoid.  The  errors  of  each  half  sinusoid  are 
uncorrolated  with  one  another,  with  amplitudes  chosen  by  a  random  number 
generator.  Figure  9  is  a  plot  of  gain  versus  phase  deviation  at  the  end  of  the 
undulator  calculated  for  different  field  errors:  (*),  no  magnetic  field  errors;  (o), 
(AB  /  B)rmg  =  0.4%;  and  (+),  (A B  /  B)rrng  =  0.5%.  The  same  seed  for  the 
random  number  generator  is  used  for  (A B / B)rm,  =  0.4%  and  (A B / B)rm,  = 
0.5%.  The  parameters  of  the  simulation  are  for  a  linearly  polarized  undulator 
with  Bu  =  5.4  kG,  A„  =  2.8  cm,  Lu  =  364  cm,  70  =  334.5,  and  e«  =  9  x  10-4  cm 
rad.  The  radiation  wavelength  is  A  =  0.25  fim.  The  beam  is  steered  on  axis  at 
the  entrance,  the  exit,  and  three  equally  spaced  axial  positions  in  the  undulator. 
The  plots  show  that  the  gain  is  inversely  related  to  the  phase  deviation.  The 
gain  degradation  is  not  significant  when  61ft  <  it/ 2. 

5.6  Error  Reduction  Techniques  and  Summary 

Several  methods  exist  for  reducing  the  detrimental  effects  of  undulator  er¬ 
rors.  In  addition  to  beam  steering,  one  may  consider  undulator  errors  that  are 
correlated.88  The  results  discussed  above  are  for  undulators  with  random  er¬ 
rors  which  are  assumed  to  be  uncorrelated  for  separation  distances  greater  than 
zac  —  Au/2.  By  considering  an  undulator  in  which  the  error  for  a  given  mag¬ 
net  pole  is  correlated  to  the  errors  of  the  surrounding  poles,  one  may  construct 
beneficial  correlations  that  reduce  the  detrimental  effects  of  the  errors. 

Alternatively,  one  may  reduce  the  detrimental  effects  of  the  errors  by  con¬ 
sidering  an  optimal  arrangement  of  the  magnet  poles.71-73  That  is,  the  magnet 
poles  are  to  be  arranged  in  such  a  way  that  the  detrimental  effects  of  the  error  of 
a  given  pole  tend  to  cancel  those  of  the  surrounding  poles.  More  specifically,  the 
magnet  poles  are  arranged  in  such  a  way  as  to  minimize  an  appropriate  “cost 
function.”  For  example,  one  may  choose  to  arrange  the  poles  such  that  the  mag¬ 
nitude  of  random  walk  \Sx\  is  minimized,  where  6x  ~  /  dz'  sin  k0(z'  -  z)6By(z'). 
(Notice  that  minimization  of  \  J  dz6B\  does  not  correspond  to  minimization  of 
|£®|.)  However,  the  results  discussed  above  indicate  that  a  more  appropriate 
cost  function  is  the  magnitude  of  the  phase  deviation  \Sift |.  Ideally,  one  would 
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like  to  maximize  the  actual  expression  for  the  gain,  Eq.  (5.6),  but  the  functional 
dependence  of  the  gain  on  the  field  errors  appears  much  too  complicated  to  be 
of  practical  use. 

The  analytical  and  numerical  work  discussed  above  indicates  that  the  phase 
deviation  Sip  is  the  single  most  important  parameter  characterizing  the  effects 
of  undulator  errors.  Although  transverse  beam  focusing  and  beam  steering  are 
highly  effective  in  controlling  the  random  walk  Sx  (in  principle,  Sx  may  be  kept 
as  small  as  desired),  this  is  not  the  case  for  the  phase  deviation  Sip.  Transverse 
beam  focusing  and  beam  steering  may  be  used  to  reduce  \Sip\  but  not  eliminate 
it.  The  phase  deviation  leads  to  a  reduction  of  FEL  gain  (the  low  gain  regime 
is  affected  more  strongly  than  the  high  gain  regime).  To  avoid  significant  loss 
in  gain,  the  above  analysis  implies  that  \Sip\  <<  27r.  Possible  error  reduction 
techniques  include  multiple  beam  steering,  correlation  of  field  errors,  and  optimal 
arrangement  of  magnet  poles.  An  optimal  arrangement  of  poles  corresponds  to 
minimization  of  |£^|,  where  Sip  is  given  by  Eq.  (5.3). 

6  OPTICAL  GUIDING 

Unlike  the  non-ideal  effects  discussed  in  Sections  3-5  that  degrade  the  per¬ 
formance  of  the  FEL,  optical  guiding  can  improve  the  gain  and  the  efficiency. 
The  FEL  mechanism  is  most  effective  when  the  radiation  beam  just  overlaps 
the  electron  beam,  such  that  the  filling  factor  /  =  oy/oy  is  approximately  unity, 
where  oy  is  the  area  of  the  electron  beam  and  oy  is  the  area  of  the  radiation. 
For  many  experiments  it  is  desirable  to  have  the  undulator  length  be  much 
longer  than  the  Rayleigh  length.  If  the  radiation  can  not  be  confined  or  guided 
by  a  waveguide  structure,  it  is  important  to  overcome  free-space  diffraction  by 
utilizing  an  important  property  of  the  FEL  -  radiation  guiding. 

In  the  one-dimensional  analysis  of  the  FEL,  the  radiation  field,  undulator 
field,  and  electron  beam  resonantly  couple  so  as  to  modify  the  longitudinal  wave 
number  of  the  radiation  field.  This  resonant  interaction,  between  the  coherent 
radiation  and  electron  beam  in  the  FEL  mechanism,  can  lead  to  focusing  of  the 
radiation  beam.  This  phenomenon  was  first  analyzed  for  the  low  gain  FEL  with 
transverse  effects,74  where  it  was  shown  that  the  diffractive  spreading  of  the 
radiation  beam  could  be  overcome  by  a  focusing  effect  arising  from  the  modified 
index  of  refraction.  Recently  optical  guiding  in  FELs  operating  in  the  expo¬ 
nential  growth  regime  has  been  studied  in  the  small  signal,  exponential  growth 
regime75-80  to  determine  the  asymptotic  behavior  of  the  radiation  beam.  Figure 
10  shows  radiation  guiding  inside  and  free-space  diffraction  outside  the  FEL  in¬ 
teraction  region  for  an  idealized  situation.  This  radiation  focusing  phenomenon 
can  play  a  central  role  in  the  practical  utilization  of  high  power  FELs.  Radiation 
self-focusing  provides  better  coupling  between  the  electron  beam  and  radiation 
beam  over  distances  longer  than  a  Rayleigh  ^ngth.  This  improves  FEL  gain 
and  efficiency  and  diminishes  losses  on  the  vacuum  chamber  walls. 
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To  illustrate  the  concept  of  radiation  self-focusing  in  the  high  gain  regime, 
consider  a  radiation  field  with  a  vector  potential  given  by 


Ar(r,2,f)  =  -A(r,z)exp[iu}(z/c  -  f)]cx/ 2  +  c.c., 


(6.1) 


where  A(r,  z)  =  \A(r,  z)\ex^r,z^  is  the  axially  symmetric  complex  amplitude  of 
the  radiation  field,  u ;  is  the  frequency,  and  c.c.  denotes  the  complex  conjugate. 
The  wave  equation  governing  Ar  is 


1  8  8  d 2  1  82  \  4tt 

t  8rV 3r  +  dz 2  c2  dt*)Ar~-  c  J*e 


(6.2) 


where  JT  is  the  driving  current  density.  Using  an  appropriate  expression  for  Jx , 
and  substituting  (6.1)  into  (6.2),  leads  to  the  following  reduced  wave  equation78 


(\  d  d 

Vrd/dr  +2t  cdz) 


a{r,z)  =  —[1  -  n2(r,2,a)]a(r,z), 


(6.3) 


where 


n(r,  z,a)  =  1  + 


1  “l{r,z) 

2  a;2 


e~*\  fBK 
To  /  \a(r,z) | 


(6.4) 


is  the  index  of  refraction  associated  with  the  medium  and  is,  in  general,  com¬ 
plex  and  a  nonlinear  function  of  a(r,  z).  Here,  u;|,(r,  z)  =  [47r|e|2n(,(r,  z)/m0]1//2 
is  the  plasma  frequency  associated  with  the  electron  beam  density,  a(r,z)  = 

\e\A/ y/2m0c2  is  the  normalized  radiation  vector  potential,  and  nb(r,z)  is  the 
electron  beam  density.  The  ()  denote  an  ensemble  average  over  all  electrons 
within  a  ponderomotive  wavelength.  The  imaginary  part  of  n  and  the  real  part 
of  n  provide  the  gain  and  the  refraction  properties  of  the  radiation,  respectively. 
Gain  and  refraction  are  coupled  in  FELs. 

Figure  11  illustrates  the  refractive  property  of  the  FEL  in  the  low  gain 
regime  when  the  input  radiation  has  a  flat  wavefront  at  the  entrance  of  the 
undulator.  Inside  the  electron  beam,  the  FEL  interaction  produces  an  index  of 
refraction  greater  than  one,  n  >  1.  The  wavefront  becomes  curved  and  radiation 
builds  up  on  the  axis. 

To  illustrate  the  property  of  radiation  self-focusing,  the  following  source- 
dependent  expansion  (SDE)78-81  for  radiation  field  a(r,z)  is  convenient: 


a{r,  z)  -  Y  aT,i (z)L, 


-[1  -  ia(r)] 


(6.5) 


where  m  =  0,1,2,  "  .  In  Eq.  (6.5),  am(z)  are  the  complex  amplitude  coefficients, 
r,(z)  is  the  radiation  spot  size,  q(z)  is  related  to  the  radius  of  curvature  of  the 
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radiation  beam  wavefront,  R  =  — (u>/2 c)r]/a  is  the  radius  of  curvature,  and  Lm 
is  the  Laguerre  polynomial. 

Let  us  assume  that  the  radiation  beam  is  approximately  Gaussian  and  the 
ao  mode  is  a  good  approximation  to  the  radiation  field.  The  unknown  quantities 
am,  r,,  and  a  can  be  solved  self-consistently.78-81 

The  beam  is  called  perfectly  guided  when  self-similar  solutions  exist,  i.e., 
the  spot  size  and  the  curvature  remain  constant.  The  complex  amplitude  of  the 
radiation  grows  exponentially, 


a(r,  z )  =  ao(z  =  0)  exp[(T  —  iAfe)z]  exp 


(6.6) 


where  T  is  the  growth  rate  and  A k  is  the  phase  shift.  Guided  radiation  beams  can 
be  obtained  in  the  small  signal,  exponential  gain  regime  for  a  uniform  electron 
beam  radius.  In  the  trapped  particle  regime,  the  radiation  spot  size  will  increase 
but  at  a  rate  slower  than  free-space  diffraction  for  the  same  beam. 

Solutions  of  r,,  a,  F,  and  A k  may  be  obtained  for  Gaussian  and  parabolic 
electron  beam  transverse  density  profiles:79-80 


n0  exp(-r2  /r2), 
no  (1  ~  r2/rl)  , 


Gaussian  electron  beam; 
parabolic  electron  beam; 


(6.7) 


where  no  is  the  electron  beam  density  on  axis.  These  solutions  are  functions  of 
the  filling  factor,  defined  as 

/  -  jb/<rr,  (6.8) 


where 


/  Gaussian  electron  beam; 

\  7T r2  /2,  parabolic  electron  beam; 


(6.9) 


and  crr  =  rrrj.  Radiation  guiding  properties  of  Gaussian  and  parabolic  electron 
beams  are  found  to  be  similar.  An  explicit  expression  for  the  filling  factor 
may  be  obtained  by  solving  an  implicit  equation  derived  in  Refs.  79  and  80. 
Consequently,  the  expressions  for  the  radius  of  curvature,  R,  the  growth  rate  T, 
the  phase  shift  Afc,  and  the  intrinsic  efficiency  i)  can  be  evaluated  directly  in 
terms  of  the  input  parameters. 

For  a  Gaussian  electron  beam  profile,  the  equation  governing  the  filling 
factor  can  be  written  as79 


where 


f,  +  f’+{\-ld)f-d2=0'  |610) 


d  = 


/  3V2kT0(rb 
^  5tt 


2/3 


(6.11) 
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To  =  (5/6fBKku)(vh)1/2(l  +  K7)~lt2 ,  (6.12) 

k  =  w/c,  and  u  =  It,(kA)/170o-  The  filling  factor  is  a  function  only  of  one 
dimensionless  parameter  d.  The  solution  for  the  filling  factor  is 


for  d  >>  1/27.  Equation  (6.13)  is  a  very  good  expression  for  the  filling  factor 
for  all  practical  parameters  of  interest. 

For  the  parabolic  electron  beam  profile,  the  equation  for  the  filling  factor 
can  be  written  as80 


2/4(4/  -  1)4(2/  -  1)  -  (6/  -  l)3d3  =  0.  (6.14) 

For  /  >  1/2,  the  approximate  solution  to  the  above  equation  is 

/  ~  ^d1//2.  (6.15) 

With  the  appropriate  explicit  expressions  for  the  filling  factor,  expressions 
can  be  determined  for  the  radius  of  curvature  R ,  the  growth  rate  I\  the  phase 
shift  Ak,  and  the  intrinsic  efficiency  q,  and  these  are  summarized  in  Table  III. 
For  large  filling  factors  (/  >  1),  the  scaling  relations  reduce  to  those  of  the  1-D 
limit.  For  /  =  1,  the  growth  rate  is  To  for  parabolic  electron  beam  profile. 

Optical  guiding  is  important  when  the  interaction  length  is  much  longer 
than  the  Rayleigh  length  and  the  Rayleigh  length  is  longer  than  the  e-folding 
length,  i.e., 

Lu  »  zr  >  1/r,  (6-16) 

where  zr  =  7rr2/A  is  the  Rayleigh  length  associated  with  the  radiation  spot  size. 
Substituting  the  expressions  from  Table  III,  we  obtain 

Ne  »  O.68(fcobr0)1/3  >  1,  (6.17) 

where  Ne  =  Lu/T  is  the  possible  number  of  e-folds  in  the  interaction  length. 
Since  it  is  desirable  to  make  the  FEL  compact,  the  FEL  undulator  length  can  be 
reduced  by  using  a  shorter  undulator  wavelength,  which  leads  to  the  requirement 
of  small  electron  beam  area  ov  To  satisfy  the  inequality  (6.17),  the  beam  current 
must  increase  as  trj,  decreases,  thus  the  beam  brightness  must  increase. 

Using  arguments  based  on  electron  trapping  in  the  ponderomotive  wave, 
the  intrinsic  efficiency  in  the  exponential  cold  Compton  regime  may  be  written 


7]  ~  Ak/ku. 


(6.18) 


The  relevant  expressions  for  efficiency  are  summarized  in  Table  III. 

The  radiation  guiding  property  of  the  FEL  can  overcome  the  reduction  of 
gain  and  efficiency  caused  by  free-space  diffraction.  For  FEL  configurations  to 
take  advantage  of  optical  guiding,  the  radiation  has  to  grow  sufficiently  fast  that 
(6.17)  is  satisfied.  At  the  same  time  high  beam  quality  is  required,  such  that  the 
effective  axial  energy  spread  is  much  less  than  the  intrinsic  efficiency,  inequality 
(3.25),  where  the  efficiency  is  defined  by  (6.18). 

7  SUMMARY 

The  above  discussion  has  delineated  the  key  physics  issues  affecting  the 
performance  of  the  FEL:  i)  electron  beam  quality,  ii)  optical  guiding  of  the 
radiation  beam,  iii)  radiation  side  band,  and  iv)  undulator  field  errors.  Perhaps 
the  most  critical  component  of  an  FEL  is  a  high  quality  electron  beam.  Several 
sources  exist  for  preducing  an  effective  energy  spread  of  the  electron  beam,  in 
addition  to  the  intrinsic  energy  spread:  transverse  beam  emittance,  undulator 
transverse  spatial  gradients,  and  undulator  field  errors.  The  electron  beam  can 
be  considered  to  be  cold  when  the  effective  axial  energy  spread  is  much  less  than 
the  intrinsic  efficiency.  Large  energy  spreads  can  significantly  reduce  the  gain 
and  efficiency.  Radiation  sidebands  and  undulator  field  errors  can  also  reduce 
the  gain  and  efficiency.  Several  methods  may  be  utilized  to  overcome  these 
detrimental  effects.  Radiation  diffraction  can  also  reduce  the  gain  and  efficiency. 
The  electron  beam,  the  radiation  field,  and  the  undulator  field  resonantly  couple 
in  the  FEL  to  modify  the  longitudinal  wave  number  of  the  radiation  field  and 
the  corresponding  index  of  refraction.  This  allows  the  radiation  beam  to  become 
“guided,”  i.e.,  the  6pot  size  r,  remains  constant  as  the  radiation  propagates. 
Radiation  guiding  provides  better  coupling  between  the  electron  beam  and  the 
radiation  beam  over  distances  longer  than  a  Rayleigh  length,  thus  improving 
both  gain  and  efficiency. 

The  physics  issues  affecting  the  performance  of  the  FEL  are  not  limited  to 
the  topics  discussed  above.  Examples  of  other  technological  and  practical  areas 
of  the  FEL  that  require  attention  are  electron  pulse  slipping  relative  to  the  ra¬ 
diation  pulse,  sensitivity  of  growth  rates  and  pulse  shape  on  mirror  detuning, 
transverse  mode  structure  and  mode  purity,  mirror  degradation  due  to  UV  ra¬ 
diation,  resonator  design  in  a  high  gain  oscillator,  electron  beam  jitter,  energy 
droop,  electron  beam  transport,  etc. 

The  FEL  is  rapidly  becoming  an  application-oriented  coherent  radiation 
source.  The  principles  of  FELs,  except  for  radiation  guiding,  have  been  demon¬ 
strated  by  the  FEL  experiments.  The  dominant  limitation  in  extending  the 
FEL  to  shorter  wavelengths  and  higher  powers  is  the  beam  quality  of  currently 
existing  accelerators.  Research  efforts  on  high  brightness  cathodes,  micro- 
undulators,  and  improved  mirrors  will  extend  the  performance  of  the  FEL  in 
the  future. 
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Table  I:  1-Dimensional  Expressions  for  Gain  and  Intrinsic  Power  Efficiency 
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Table  II:  Accelerators  Used  in  FELs3'5’82  88 
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9  This  is  the  core  brightness.  The  edge  brightness  would  be  approximately  a  factor  of  2  lower. 
*  Measured  normalized  rms  emittance  is  Tren<rrri,  =  6  x  10-2  cm  rad.  Used  tn  =  3en>rm, 
and  Bn  -  2Ih/9n2e2n  PmJ. 

’  Same  for  ail  energies. 

;  Design  parameters. 


Table  III:  Explicit  Conditions  for  Guided  Radiation  Beam  in  the 
Compton  Exponential  Gain  Regime 
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Budker’s  parameter 
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electron  beam  radius 
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Reflector 
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Fig.  1  —  Schematic  of  an  FEL  consisting  of  an  electron  beam,  an  undulator  and  the  radiation,  in  an 
oscillator  configuration 
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Distribution  function  of  electrons  at  the  entrance  and  exit  of  a  tapered  undulator 
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z-°  z-Lu 

(b) 


Fig.  7  —  (a)  Radiation  amplitude  modulated  by  sidebands,  and  (b)  oscillation  of  the  separatrix  as  a 
result  of  the  modulated  radiation  amplitude 


41 


NORMALIZED  GAIN,<G> 


FREQUENCY  MISMATCH,  fi 

Fig.  8  —  Mean  normalized  gain,  <G/,  versus  frequency  mismatch,  n.  for  several  values  of  rms  field 
error,  6 for  a  linearly  polarized  undulator  with  Bu  =  5.4  kG,  Xu  =  2.8  cm.  Lu  =  3.6  m,  and 
7  =  350  in  the  limit  kB  =  0 


42 


1 


1 


^  1 

s 

1 

c 

’o 

o 


<D 

$ 

o 

CL 


0  1.0  2.0  3.0  4.0  5.0 

Phase  Deviation  <5^( z  =  Lu) 

Fig.  9  —  Plot  of  gain  versus  phase  deviation  for  different  field  errors:  (*)  no  magnetic  field  errors, 
(:)  (A B/B)^  =  0.4%,  and  (  +  )  ( AB /fl )rms  =  0.5%.  The  same  random  number  seeds  are  used 
for  ( : )  and  ( + ) 
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— ]  FEL  Interaction  Region  p— 

Free  Space*  *  Free  Space 

Diffraction  Optical  Guiding  j  Diffraction 

Fig.  10  —  Schematic  of  an  FEL  where  the  radiation  is  “perfectly”  guided  in  the  interaction  region 
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illustrating  the  refractive  property  of  the  FEL  in  the  low  gain  regime 
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